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Abstract 

Charmonium production and suppression in In+In and Pb+Pb reactions at SPS 
energies is investigated with the HSD transport approach within the 'hadronic co- 
mover model' as well as the 'QGP threshold scenario'. The results of the transport 
calculations for J/^> suppression and the ^f' to J/^> ratio are compared with the 
recent data of the NA50 and NA60 Collaborations. We find that the comover ab- 
sorption model - with a single parameter |Mo| 2 for the matrix element squared for 
charmonium-meson dissociation - performs best with respect to all data sets. The 
'threshold scenario' - within different assumptions for the melting energy densities 
- yields a reasonable suppression for J /iff but fails in reproducing the iff' to J /\P 
ratio for Pb+Pb at 158 A-GeV. Predictions for Au+Au reactions are presented for 
a bombarding energy of 25 A-GeV in the different scenarios which will allow for a 
clear distinction between the models from the experimental side at the future FAIR 
facility. 
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1 Introduction 



The dynamics of ultra-relativistic nucleus-nucleus collisions at Super-Proton- 
Synchrotron (SPS) and Relativistic-Heavy-Ion-Collider (RHIC) energies are 
of fundamental interest with respect to the properties of hadronic/partonic 
systems at high energy densities. Especially the formation of a quark-gluon 
plasma (QGP) and its transition to interacting hadronic matter has motivated 
a large community for almost three decades [T] . The c, c quark degrees of 
freedom are of particular interest in context to the phase transition to the 
QGP, since cc meson states might no longer be formed due to color screening 
[2T314T5] . However, more recent lattice QCD (1QCD) calculations have shown 
that the J/\I/ survives up to at least 1.5 T c (T c ^170 MeV) such that the 
lowest cc states remain bound up to rather high energy density |6,7 8j. On 
the other hand the \c and appear to melt soon above T c . It is presently 
not clear if also the D or D* mesons will survive at temperatures T > T c but 
strong correlations between a light quark (antiquark) and a charm antiquark 
(quark) are likely to persist [9]. 

The standard approach to charmonium production in heavy-ion collisions as- 
sumes that cc pairs are created exclusively at the initial stage of the reaction in 
primary nucleon-nucleon collisions. At the very early stage color dipole states 
are expected to be formed which experience i) absorption by interactions with 
nucleons of the colliding nuclei (c/. Refs. [TUfTI] ). These cc states are assumed 
to be absorbed in a 'pre-resonance state' before the final hidden charm mesons 
are formed. This absorption - denoted by 'normal nuclear suppression' - is 
also present in p+A reactions and determined by a dissociation cross section 
Ob ~ 4 to 7 mb). Those charmonia or 'pre-resonance' states - that survive 
normal nuclear suppression - furthermore suffer from ii) a possible dissocia- 
tion in the deconfined medium at sufficiently high energy density and iii) the 
interactions with secondary hadrons (comovers) formed in a later stage of the 
nucleus-nucleus collision. 

The geometrical Glauber model of Blaizot et al. [12] , as well as the percolation 
model of Satz |4] , assumes that the QGP suppression ii) sets in rather abruptly 
as soon as the energy density exceeds a threshold value e c , which is a free 
parameter. This version of the standard approach will be referred to as the 
QGP 'threshold scenario'. The latter model is motivated by the idea that the 
charmonium dissociation rate is drastically larger in a quark-gluon-plasma 
(QGP) than in a hadronic medium |!4J such that further comover absorption 
channels might be neglected. 

On the other hand, the extra suppression of charmonia in the high density 
phase of nucleus-nucleus collisions at SPS energies |13|I14II15II16|I17II18| has been 
attributed to inelastic comover scattering (c/. |Hlll9f2Uf21ll22f23f24ll25] and 
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Refs. therein) assuming that the corresponding J/\I/-hadron cross sections are 
in the order of a few mb [26 27.28 ,29j . Theoretical estimates here differ by 
more than an order of magnitude [30j especially with respect to J/\l/-meson 
scattering such that the question of charmonium suppression is still open. 
Additionally, alternative absorption mechanisms - such as gluon scattering on 
color dipole states - might play a role as suggested in Refs. [3Tll32f33|l3"4"j and 
also lead to a reduction of the final J/\I/ formation in central nucleus- nucleus 
collisions. 



We recall that apart from absorption or dissociation channels for charmonia 
also recombination channels such D + D — > J/ty + meson may play a role. A 
previous analysis within the HSD transport approach [35J - employing the co- 
mover absorption model - has demonstrated that the charmonium production 
from open charm and anticharm mesons indeed becomes essential in central 
Au+Au collisions at RHIC. This is in accordance with independent studies in 
Refs. [281132] and also with the data from PHENIX [36]. On the other hand, 
these backward channels - relative to charmonium dissociation with comoving 
mesons - have been found to be practically negligible at the SPS energies of 
interest here. Nevertheless, in our dynamical studies below we will include the 
'backward' channels for completeness. 



A couple of models have predicted J/\l/ suppression in In+In collisions as a 
function of centrality at 158 A-GeV based on the parameters fixed for Pb+Pb 
reactions at the same bombarding energy. However, the predictions within 
the comover model and 'threshold scenario' from Refs. [38I37|39] have failed 
to describe the data with sufficient accuracy. This might be either due to 
the missing dynamics of the nucleus-nucleus collisions in these models or to 
inadequate physical assumptions about the dissociation mechanism. 



In the present work we extend the previous studies within the comover model 
in Refs. [4T)|35|I41] and test the 'QGP threshold scenario' - described in Section 
2 - in comparison to the Pb+Pb data at 158 A-GeV from NA50 as well the 
high statistics data from NA60 for In+In collisions at the same bombarding 
energy. The question we aim at solving in Section 3 is: 1) can any of the 
models be ruled out by the combined data sets and 2) do the recent NA60 
data provide a hint to QGP formation at the top SPS energy? In Section 
4 we, furthermore, will provide predictions for the charmonium suppression 
in Au+Au collisions at 25 A-GeV that will be measured at the future FAIR 
facility by the CBM Collaboration. 
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2 Brief description of charmonium channels in HSD 

The microscopic Hadron- String-Dynamics (HSD) transport calculations (em- 
ployed here) provide the space-time geometry of the nucleus-nucleus reaction 
and a rather reliable estimate for the local energy densities achieved since the 
production of secondary particles is described rather well from SIS to RHIC en- 
ergies [12] . In order to examine the dynamics of open charm and charmonium 
degrees of freedom during the formation and expansion phase of the highly 
excited system created in a relativistic nucleus-nucleus collision within trans- 
port approaches, one has to know the number of initially produced particles 
with core quarks, i.e. D, D, D\ D\ D s , D s , D*, D* s , J/V(1S),V'(2S), Xc(lP)- 
In this work we follow the previous studies in Refs. [T9|22|I35||4"0] and fit the 
total charmonium cross sections (i = Xc, J ,^') from NN collisions as a 
function of the invariant energy ^fs by the expression 

•rM = /.«(i-^)°OV-ViS), (i) 

where denotes the mass of charmonium i while ^/soi = rrii + 2rriN is the 
threshold in vacuum. The parameters in ([1]) have been fixed to describe the 
J/^ and ^' data up to the RHIC energy y/s = 200 GeV (cf. Fig. □). We 
use a = 0.2 mb, a — 10, /3 = 0.775. The parameters fi are fixed as f Xc = 
0.636, fj/i& = 0.581, f^i = 0.21 in order to reproduce the experimental ratio 

gto - J/9)c M +_gfag -~ J/»K. = .3 44 ± o. 031 

measured in pp and ttN reactions |4*3|4*4] as well as the averaged pp and pA 
ratio (B ^{^>')a^i) / (B ^(J /^)a j ~ 0.0165 (cf. the compilation of experi- 
mental data in Ref. [45J). Here the experimentally measured J/\l/ cross section 
includes the direct J/\l/ component (crj/^) as well as the decays of higher char- 
monium states Xci ^' i i- e - 

v e J/l = °JI* + B( X c - J/*)<r Xc + B(V - J/*)<T 9 ,. (2) 

Note, we do not distinguish the Xci(l-P) and Xc2(lP) states. Instead, we use 
only the Xci(l-P) state (which we denote as Xc), however, with an increased 
branching ratio for the decay to J/^f in order to include the contribution of 
X c2 (lP), ie. B( X c -> J/V) = 0.54. We adopt B(V -> J/V) = 0.557 from 
Ref. 

In addition to primary hard NN collisions the open charm mesons or charmo- 
nia may also be generated by secondary meson-baryon (mB) reactions. Here 
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Fig. 1. The cross section for D + D, J/\I/ and VP' meson production in pN 
(left part) and irN reactions (right part). The solid lines show the parametri- 
sations used in HSD, whereas the symbols stand for the experimental data 
|48|49|50|51|52|53|54|55|56|57j . The J/^f cross sections include the decay from Xc 
mesons. 

we include all secondary collisions of mesons with baryons by assuming that 
the open charm cross section (from Section 2 of Ref. [ID]) only depends on 
the invariant energy \fs and not on the explicit meson or baryon state. Fur- 
thermore, we take into account all interactions of 'formed' mesons - after a 
formation time of Tp = 0.8 fm/c (in their rest frame) j3?J ~~ with baryons or 
diquarks, respectively. For the total charmonium cross sections from meson- 
baryon (or ttN) reactions we use the parametrization (in line with Ref. |20j): 



*I N (s) = fi b [ 1 - ^ ] (3) 




with 7 = 7.3 and b = 1.24 mb, which describes the existing experimental data 
at low y/s reasonably well as seen from Fig. HJ 

Apart from the total cross sections, we also need the differential distribution 
of the produced mesons in the transverse momentum pp and the rapidity y (or 
Feynman Xp) from each individual collision. We recall that xp = p z /p™ ax ~ 
2p z /\/s with p z denoting the longitudinal momentum. For the differential 
distribution in xp from iViV and ttN collisions we use the ansatz from the 
E672/E706 Collaboration [58]: 

dN 

> (1 - \x F \) c exp(-b PT p T ), (4) 



dxpdp' 



where b PT = 2.08 GeV -1 and c = a/(l + b/ y/s). The parameters a, b are chosen 
as a NN = 13.5, bj^N = 24.9 for iVTV collisions and a nN = 4.11, b wN = 10.2 for 
ttN collisions as in [35~40j. 

The parametrizations of the total and differential cross sections for open charm 
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mesons from pN and nN collisions are taken as in Refs. (3500]. Here we show 
only the total cross sections for D + D productions in Fig. [TJ 

In order to study the effect of charmonium rescattering we adopt the follow- 
ing dissociation cross sections of charmonia with baryons independent of the 
energy (in line with the recent NA50 and NA60 compilations [TBlloT?] ): 

a c - cB = 4.18 mb; (5) 
o~j/^b = 4.18 mb; <j Xc b = 4.18 mb; ays = 7.6 mb. 

In (jSJ) the cross section a c5 B stands for a (color dipole) pre-resonance (cc) - 
baryon cross section, since the cc pair produced initially cannot be identified 
with a particular hadron due to the uncertainty relation in energy and time. 
For the life-time of the pre-resonance cc pair (in it's rest frame) a value of r cg 
= 0.3 fm/c is assumed following Ref. [60]. This value corresponds to the mass 
difference of the 1 J r ' and J/^>. 

For D, D*, D, D* - meson (ir, rj, p, u) scattering we address to the calculations 
from Ref. [271128] which predict elastic cross sections in the range of 10-20 mb 
depending on the size of the formfactor employed. As a guideline we use a 
constant cross section of 10 mb for elastic scattering with mesons and also 
baryons, although the latter might be even higher for very low relative mo- 
menta. Since the .D-meson dynamics is of minor importance for charmonium 
regeneration at SPS energies we discard a more detailed description. 



2.1 The comover absorption model 



As already pointed out before, the J/^f, Xc formation cross sections by open 
charm mesons or the inverse comover dissociation cross sections are not well 
known and the significance of these channels is discussed controversely in the 
literature [30|62|l63|l64ll65|l66|l67j . We here follow the concept of Refs. [3"oT4"Tj and 
introduce a simple 2-body transition model with a single parameter |M | 2 , that 
allows to implement the backward reactions uniquely by employing detailed 
balance for each individual channel. 

Since the charmonium-meson dissociation and backward reactions typically 
occur with low relative momenta ('comovers') it is legitimate to write the 
cross section for the process 1 + 2^3 + 4 as 

/ x i E i E 2E 3 E A ~ 2 fm 3 + mA 6 p f 
^1+2^3+4(3) = 2 — -i-, (6) 

8 \ V s J Pi 



where E k denotes the energy of hadron k (k — 1,2,3,4), respectively. The 



6 



initial and final momenta for fixed invariant energy y^s are given by 



2 _ (s - (mi + m 2 ) )(s - (mi - m 2 ) ) 
Pi ~ Is ' 

p2 = (g - ( m 3 + m 4 ) 2 )(s - (m 3 - m 4 ) 2 ) ^ 

where m^ denotes the mass of hadron k. In ([6]) |Mj| 2 (i = Xc J/i>, stands for 
the effective matrix element squared, which for the different 2-body channels 
is taken of the form 



\Mi\ 2 = \Mi\ 2 for (it, p) + (cc)i —> D + D (8) 
\M i \ 2 = 3\M i \ 2 for (tt, p) + (cc), ^ D* + D, D + D*, D* + D* 

\M l \ 2 = ^\M l \ 2 for (K,K*) + (cc) t ^ D s + D, D s + D 

\Mi\ 2 = \Mi\ 2 for (K,K*) + (cc)i -> D s + D* , D s + D*, D* s +D, 

D* s +D, D* s +D* 

The relative factors of 3 in (jSJ) are guided by the sum rule studies in [68] which 
suggest that the cross section is increased whenever a vector meson D* or D* 
appears in the final channel while another factor of 1/3 is introduced for each 
s or s quark involved. The factor ((m.3 + m^/y/s) 6 in ([6]) accounts for the 
suppression of binary channels with increasing yfs and has been fitted to the 
experimental data for the reactions tc + N p + N,u + N,<f) + N, K + + A in 
Ref. |69]. 

We use (for simplicity) the same matrix elements for the dissociation of all 
charmonium states i (i = Xc J/' l P,' l l J ') with mesons: 



\Mj,^\ 2 = \M x f = |M*f = |M | 2 . (9) 

We note that in Ref. [33] the parameter |M | 2 was fixed by comparison to 
the J/ 1 ® suppression data from the NA38 and NA50 Collaborations for S+U 
and Pb+Pb collisions at 200 and 158 A-GeV, respectively. In the present 
study however, this parameter has to be readjusted in accordance with the 
updated value of the cross section <^ of charmonium dissociation on baryons 
(following the latest NA50 and NA60 analysis [T5|5^] ). The best fit is obtained 
for |M | 2 = 0.18 fm 2 /GeV 2 . 

The advantage of the model introduced in [35fH] is that detailed balance for 
the binary reactions can be employed strictly for each individual channel, i.e. 
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(2S 1 + 1)(2S 2 + 1) pf 
(25 , 3 + l)(25 , 4+l) p}' 



^3+4^1+2(s) = ^l+2^3+4(g) /00 , — — , (10) 



and the role of the backward reactions ((cc)j+meson formation by £) + D 
flavor exchange) can be explored without introducing any additional param- 
eter once \Mi\ 2 is fixed. In Eq. ( TTDT) the quantities Sj denote the spins of the 
particles, while pf and p 2 denote the cms momentum squared in the initial 
and final channels, respectively. The uncertainty in the cross sections (11 01) 
is of the same order of magnitude as that in Lagrangian approaches using 
e.g. SU (4) fi avor symmetry [2Tf2g] . since the formfactors at the vertices are 
essentially unknown [BE]- It should be pointed out that the comover dissocia- 
tion channels for charmonia are described in HSD with the proper individual 
thresholds for each channel in contrast to the more schematic comover absorp- 
tion model [TT] . 

We recall that (as in Refs. [35|40f70|l71f72] ) the charm degrees of freedom in 
the HSD approach are treated perturbatively and that initial hard processes 
(such as cc or Drell-Yan production from NN collisions) are 'precalculated' to 
achieve a scaling of the inclusive cross section with the number of projectile 
and target nucleons as Ap x Ax when integrating over impact parameter. For 
fixed impact parameter h the cc yield then scales with the number of binary 
hard collisions N bin (cf. Fig. 8 in Ref. [3D]). 



2.2 Implementation of the 'threshold scenario' 



The HSD transport model allows to calculate the energy-momentum tensor 
T flv (x) for all space-time points x and thus the energy density e(x) in the local 
rest frame. In order to exclude contributions to from noninter acting nu- 
cleons in the intial phase all nucleons without prior interactions are discarded 
in the rapidity intervals [y tar — 0.4, y tar + 0.4] and [y pro — 0.4, y pro + 0.4] where 
y tar and y pro denote projectile and target rapidity, respectively. Note that the 
initial rapidity distributions of projectile and target nucleons are smeared out 
by about ±0.4 due to Fermi motion. 

In the actual calculation the initial grid has a dimension of 1 fm x 1 fm 
x l/jcm fm, where 7 cm denotes the Lorentz 7-factor in the nucleon-nucleon 
center-of-mass system. After the time of maximum overlap t m of the nuclei 
the grid-size in beam direction Az = l/7 cm [fm] is increased linearly in time 
as Az = Az + a(t — t m ), where the parameter a is chosen in a way to keep 
the particle number in the local cells roughly constant during the longitudinal 
expansion of the system. In this way local fluctuations of the energy density 
e(x) due to fluctuations in the particle number are kept low. 

As an example we display in Fig. [2] the energy density e(x, y = 0, z;t) for a 
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e(x,y=0,z;t) [GeV/fm 3 ] 



Pb+Pb @ 160 A GeV, b=l fm 




-4 -2 2 4 -4 -2 2 4 




^4 ' ^2 ' ' 2 4 -4 ' ~-2 ' ' 2 ' 4 

z [fm] z [fm] 



Fig. 2. The energy density e(x, y = 0, z; i) from HSD for a Pb+Pb collision at 160 
A-GeV and impact parameter b = 1 fm in terms of contour lines (0.01, 1, 2, 3, 4 
GeV/fm 3 ) for times of 1, 2, 3 and 5 fm/c (from contact). Note that noninteracting 
nucleons have been discarded in the actual calculation of the energy-momentum 
tensor. 

Pb+Pb collision at 160 A- GeV and impact parameter b = 1 fm in terms of 
contour lines for times of 1, 2, 3 and 5 fm/c (from contact). It is clearly seen 
that energy densities above 4 GeV/fm 3 are reached in the early overlap phase 
of the reaction and that e(x) drops within a few fm/c below 1 GeV/fm 3 in 
the center of the grid. On the other hand the energy density in the region of 
the leading particles - moving almost with the velocity of light - stays above 1 
GeV/fm 3 due to Lorentz time dilatation since the time t here is measured in 
the nucleon-nucleon center-of-mass system. Note that in the local rest frame of 
the leading particles the eigentime r is roughly given r « t/j cm with 7 cm ^ 9.3. 

Another view of the space time evolution of the energy density is given in 
Fig. [3] where we display e(x — 0, y = 0, z\ t) for the same system as in Fig. [2] 
on a linear scale. The contact time of the two Pb nuclei here is 2 fm/c and the 
overlap phase of the Lorentz contracted nuclei is identified by a sharp peak in 
space-time which is essentially given by the diameter of the nuclei divided by 
7 cm . As noted before, the energy density in the center of the reaction volume 
(z ~ 0) drops fast below 1 GeV/fm 3 whereas the ridges close to the lightcone 
basically stem from the leading ends of the strings formed in the early nucleon- 
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Fig. 3. The energy density e(x = 0, y = 0, z; t) from HSD for a Pb+Pb collision at 
160 A-GeV and impact parameter b = 1 fm on a linear scale. Note that noninteract- 
ing nucleons have been discarded in the actual calculation of the energy-momentum 
tensor such that e(x) 7^ only after contact of the two Pb nuclei which is ~ 2 fm/c. 

nucleon collisions. In these space-time regions all reaction rates are reduced 
by the factor ~ 1/jcm such that the transport calculations have to be carried 
to large times of several hundred fm/c in order to catch the dynamics and 
decays in these regions. In the central regime, however, all interaction rates 
vanish after about 15 fm/c. Since the c, c pairs are produced dominantly at 
midrapidity with a small spread in rapidity {a y ~ 0.8 at 160 A-GeV) it is the 
central region that is of primary interest for this study. 

Of further interest is the size of the total volume (measured in the nucleon- 
nucleon cms) with an energy density above a certain cut e c , i.e. 

V(e c ;t):= J d 3 r 0(e(r; t) - e c ), (11) 



which quantifies the volume for charmonium dissolution as a function of time t. 
The corresponding information is displayed in Fig. H]for e c = 1 GeV/fm 3 (left 
part) and 1.5 GeV/fm 3 (right part) as a function of time for impact parameter 
b=l to 12 fm (in steps of A b = 1 fm). These volumina may be compared 
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i ^ Pb+Pb, 160 A GeV | 



H Pb+Pb, 160 A GeV | - 



e > 1.5 GeV/fm 
b=l,2,3, 12 fm 




5 10 

time [fm/c] 



5 10 

time [fm/c] 



Fig. 4. The volume V{e c ;t) (JTTJ) from HSD for Pb+Pb collisions at 160 A -GeV and 
impact parameter b = 1,2, .., 12 fm for e > e c = 1 GeV/fm 3 (left part) and e > e c 
= 1.5 GeV/fm 3 (right part). 

to the Lorentz contracted eigenvolume of a Pb nucleus that is about 160 fm 3 
in the cms. It is clearly seen that hadron formation and the explosion (or 
expansion) of the system lead to larger volumina V(e c ; t) in central reactions 
especially for e c = 1 GeV/fm 3 . Note, however, that charmonia dynamically 
cannot explore the whole volume displayed in Fig. HI since this volume is 
dominated by the space-time regimes close to the lightcone (cf. Fig. [3]), where 
practically no charmonia appear at 160 A-GeV. 

The 'threshold scenario' for charmonium dissociation now is implemented in 
a straight forward way: whenever the local energy density e(x) is above a 
threshold value Ej, where the index j stands for J/^f, Xc, the charmonium 
is fully dissociated to c + c. The default threshold energy densities adopted are 
e 1 = 16 GeV/fm 3 for J/tf , e 2 = 2 GeV/fm 3 for Xc, and e 3 = 2 GeV/fm 3 for 
The reformation of charmonia at the phase boundary to the hadronic system 
is discarded in view of the very low charm quark density at SPS energies. 



3 Comparison to data 



We directly step on with results for the charmonium suppression at SPS en- 
ergies in comparison with the experimental data from the NA50 and NA60 
Collaborations. These Collaborations present their results on J/ty suppression 
as the ratio of the dimuon decay of J/^> relative to the Drell-Yan background 
from 2.9 - 4.5 GeV invariant mass as a function of the transverse energy Et-i 
or alternative, as a function of the number of participants A^art, i.e. 

B^a(J/*)/a(DY)\ 2 , 9 ^. 5 , (12) 
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Fig. 5. The ratio B^^a(J/^f)/a(DY) as a function of the number of participants in 
In+In (l.h.s.) and Pb+Pb reactions (r.h.s.) at 158 A-GeV. The full symbols denote 
the data from the NA50 and NA60 Collaborations (from Refs. |17ll8|61j ), while 
the dashed (blue) lines represent the HSD calculations including only dissociation 
channels with nucleons. The lower parts of the figure show the HSD results in the 
same limit for the 'J?' to J / ^ ratio as a function of N par t (for In+In) or the transverse 
energy Et (for Pb+Pb). The solid (red) lines show the HSD results for the comover 
absorption model with a matrix element squared |Mo| 2 = 0.18 fm 2 /GeV 2 . The (light 
blue) bands in the upper parts of the figure give the estimate for the normal nuclear 
J/'i? absorption as calculated by the NA60 Collaboration. The vertical lines on the 
graphs reflect the theoretical uncertainty due to limited statistics of the calculations. 



where 5 W is the branching ratio for J/ty — > . 

In the theoretical approaches we calculate the J/^/ survival probability Sj/q, 
denned as 

N fln 

s J/f - —jyyi ( 13 ) 
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Fig. 6. Same as Fig.[5]but for the 'QGP threshold scenario' with £j/^ = 16 GeV/fm 3 , 
e Xc = 2 GeV/fm 3 = e^i while discarding comover absorption, i.e. for |Mo| 2 = 0. 

where N J J^ and N J Jg denote the final number of J/\P mesons and the num- 
ber of J/^'s produced initially by BB reactions, respectively. In order to 
compare our calculated results to experimental data we need an extra in- 
put, i.e. the normalization factor B^unn^J 7^) / 'ctnn(DY), which defines the 
J/ty over Drell-Yan ratio for elementary nucleon-nucleon collisions. We choose 
B^<J NN (J/^)/<T NN (DY) = 36 in line with the NA60 compilation [T8] . 



Furthermore, the suppression is presented experimentally by the ratio 

B^> -> nn)a{W)/a{DY) 
B^(J/V - w)<t(J/*)/<t(DY)' 



(14) 



In our calculations we adopt this ratio to be 0.0165 for nucleon-nucleon colli- 
sions, which is again based on the average over pp,pd,pA reactions |45j. 

In order to investigate the anomalous charmonium suppression in nucleus- 
nucleus collisions we first show in Fig.[5]the calculated ratio B flll a(J/^f)/a(DY) 
as a function of N pait for Pb+Pb and In+In collisions at 158 A- GeV (upper 
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plots) in the nuclear suppression scenario, i.e. without comover dissociation 
or 'QGP threshold suppression'. The dashed (blue) lines stand for the HSD 
result while the (light blue) bands give the estimate for the normal nuclear 
J/ty absorption as calculated by the NA60 Collaboration. The normal nuclear 
suppression from HSD is seen to be slightly lower than the (model dependent) 
estimate from NA60, however, agrees quite well with their model calculations 
for more central reactions. The various experimental data points have been 
taken from Refs. [T7|18|61j . It is clearly seen that the charmonium dissociation 
with only nucleons is insufficient to describe the data for both systems. In the 
lower part of Fig. Owe compare the calculated ratio ^ over J/ty as a func- 
tion of Np ar t (for In+In reactions) or the transverse energy Et (for Pb+Pb 
collisions), respectively, in comparison to the data available. The Pb+Pb data 
demonstrate that the centrality dependence as well as the absolute ratio can- 
not be explained by nuclear dissociation channels alone which is a well known 
fact in the community. 

Apart from the statistical uncertainties in the calculations - reflected by the 
vertical lines on the theoretical graphs in Fig. [5] - some dependence on the 
model parameters enters the actual numbers in Fig. The charmonium nu- 
clear absorption cross section is considered to be 'fixed' by the NA50/NA60 
compilations and we have taken the same cross section for the 'pre-resonance' 
cross section for the J/^f and Xc- Accordingly the life-time of the pre-resonance 
state (r C c = 0.3 fm/c) has no impact on the absorption with baryons as far 
as the J/ty and Xc mesons are concerned. Only for collisions with baryons 
this plays a role since the + baryon cross section is larger (7.6 mb). Con- 
sequently the J/ty suppression (including the feed down from Xc) does not 
depend on r cE . The finite life-time, however, plays a role for suppression as 
can be seen in the lower right part of Fig. [5] since it leads to a larger baryon 
absorption of (relative to J/^) with increasing centrality by about 22%. 
According to our understanding the life-time of 0.3 fm/c is a lower limit (in 
line with the Heisenberg uncertainty relation); it leads to a maximum sup- 
pression of ty' relative to J/^f with centrality. On the other hand, a very large 
life-time t c5 will lead to a constant ratio of to J/ty with centrality since only 
the pre-resonance cross section will apply. Accordingly, the to J/ty ratio is 
driven by the value of t c5 since the ratio of the dissociation cross section for 
the formed mesons is fixed by the ratio of their mean square radii. However, 
independently on the life-time r ce of the pre-resonance state the experimental 
data will be badly missed for the to J/^f ratio if only baryon dissociation 
is included, because the case considered here already provides a maximum 
suppression of the for interactions with baryons. 

As a next step we add the comover dissociation channels within the model 
described in Section 2.1 for a matrix element squared |Mo| 2 = 0.18 fm 2 /GeV 2 . 
Note that in this case the charmonium reformation channels are incorporated, 
too, but could be discarded since the charmonium regeneration is negligible at 
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SPS energies (cf. Ref. [35]). The extra suppression of charmonia by comovers 
is seen in Fig. [5] (solid (red) lines) to match the J/ty suppression in In+In 
and Pb+Pb as well as the ty' to J/\I/ ratio (for Pb+Pb) rather well. The more 
recent data (1998-2000) for the ty' to J/^> ratio agree with the HSD prediction 
within error bars. This had been a problem in the past when comparing to the 
1997 data (dark green stars). We conclude that the comover absorption model 
presently cannot be ruled out on the basis of the available data sets within 
error bars. The \&' to J/^> ratio for In+In versus centrality is not yet available 
from the experimental side but the theoretical predictions are provided in Fig. 
Eland might be approved/falsified in near future. 

Some comments to the comover absorption model appear in place: As shown 
in Fig. 7.2 of Ref. [19] the comover densities in central Pb+Pb collisions at 
158 A-GeV become quite large and almost reach 2/fm 3 in the maximum which 
appears high for 'free' mesons with an eigenvolume of about 1 fm 3 . However, 
the quasi-particle mesons considered here dynamically should not be identi- 
fied with 'free' meson states that show a long polarization tail in the vacuum. 
As known from lattice QCD the correlators for pions and p-mesons survive 
well above the critical temperature T c , such that 'dressed' mesons, i.e spectral 
densities with the quantum numbers of the pseudo-scalar and vector (isovec- 
tor) modes, also show up at high energy density (similar to the J/^f discussed 
above [6l7,8j). Such 'dressed' mesons are expected to have a shorter polariza- 
tion tail - since the reference vacuum has changed and the vacuum polarization 
decreases - and thus are much smaller in size. In any case, these 'states' (or 
resonances) will dissociate charmonia by a quark rearrangement interaction in 
exchanging a light quark with a c or c quark. 

The results for the 'threshold scenario' are displayed in Fig. [6] in comparison 
to the same data for the thresholds ej/q, = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 
= e^i while discarding any dissociation with comovers, i.e. |M | 2 =0. In this 
scenario the J/^> suppression is well described for In+In but the suppression is 
slightly too weak for very central Pb+Pb reactions. This result emerges since 
practically all Xc and dissolve for N part > 100 in both systems whereas the 
J/ty itself survives at the energy densities reached in the collision. Since the 
nucleon dissociation is a flat function of N part for central reactions the total 
absorption strength is flat, too. The deviations seen in Fig. [6] might indicate 
a partial melting of the J/\I/ for N part > 250, which is not in line with the 
lattice QCD calculations claiming at least Ej/q, > 5 GeV/fm 3 . In fact, a lower 
threshold of 5 GeV/fm 3 (instead of 16 GeV/fm 3 ) for the J/^> has practically no 
effect on the results shown in Fig. [6J Furthermore, a threshold energy density 
of 2 GeV/fm 3 for the leads to a dramatic reduction of the to J/^> ratio 
which is in severe conflict with the data (lower part of Fig. [6]). Also note that 
due to energy density fluctuations in reactions with fixed N part (or Er) there 
is no step in the suppression of J/ ^ versus centrality as pointed out before by 
Gorenstein et al. in Ref. [75] . 
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Fig. 7. Same as Fig.[5]but for the 'QGP threshold scenario' with ej/q = 16 GeV/fm 3 , 
e Xc = 2 GeV/fm 3 , e^i = 6.55 GeV/fm 3 while discarding comover absorption, i.e. 
for |M | 2 = 0. 

Since in the 'threshold scenario' the J/\l/ suppression is rather well accounted 
for by the melting of the Xc it might be tempting to 'extract' a dissociation 
energy density for the via the to J/^> ratio. This may indeed be achieved 
for e^i = 6.55 GeV/fm 3 as shown in Fig. [7] where now also the to J/\l/ ratio 
is well reproduced for Pb+Pb. Respective predictions for In+In within this 
scenario are presented in the lower left part of Fig. [7| and wait for confirmation 
or disproof. We recall, however, that the energy density of 6.55 GeV/fm 3 for 
ty' dissociation is not supported by present lattice calculations such that this 
limit should be ruled out. 

As the last model scenario we assume that both J/\l/ and are practically 
not dissociated in nucleus-nucleus reactions at SPS energies but let the Xc 
melt above 2 GeV/fm 3 as before. In this case one may estimate the maximum 
suppression by comovers that is compatible with the data. This combined 
scenario allows to fix a matrix element squared |M | 2 ~ 0.09 fm 2 /GeV 2 that 
gives results still compatible with the J/\l/ suppression for In+In and Pb+Pb 
(solid lines in Fig. El). The dashed lines in Fig. M show the results without any 
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Fig. 8. Same as Fig.[5]but for the 'QGP threshold scenario' with Ejm = 16 GeV/fm 3 , 
e Xc = 2 GeV/fm 3 , = 16 GeV/fm 3 (dashed blue lines) and with additional 
comover absorption for |Mo| 2 = 0.09 fm 2 /GeV 2 (red solid lines). 

comover absorption and are due to the melting of the Xc above 2 GeV/fm 3 
alone. However, in this combined approach the to J/\l/ ratio turns out to be 
systematically too high in comparison to the data for Pb+Pb such that this 
scenario should be ruled out, too, in particular with respect to the extreme 
threshold for ty' dissociation. 

Since the NA60 Collaboration prefers to represent their data in a model depen- 
dent way by plotting their experimental results relative to the normal nuclear 
absorption model we additionally show in Fig. [9] our calculations for In+In 
(red lines with open squares) and Pb+Pb (blue lines with open circles) as a 
function of the number of participants N part relative to the normal nuclear 
absorption given by the straight black line (according to the NA60 compila- 
tion). The full dots and squares denote the respective data from the NA50 
and NA60 Collaborations. The model calculations reflect the comover absorp- 
tion model (right part) and the 'QGP threshold scenario' (left part) with 
Ej^ = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 , e«* = 6.55 GeV/fm 3 . Since only the 
representation is different the message stays the same: The comover absorp- 
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Fig. 9. The ratio B flfl a(J/^)/a(DY) as a function of the number of participants 
N par t in In+In (red line with open squares) and Pb+Pb reactions (blue line with 
open circles) at 158 A-GeV relative to the normal nuclear absorption given by 
the straight black line. The full dots and squares denote the respective data from 
the NA50 and NA60 Collaborations. The model calculations reflect the comover 
absorption model (right part) and the 'QGP threshold scenario' (left part) with 
ejjq, = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 , eyi = 6.55 GeV/fm 3 while discarding 
comover absorption. 

tion model follows slightly better the fall of the J/^f survival probability with 
increasing centrality whereas the 'threshold scenario' leads to an approximate 
plateau in both reactions for high centrality. 



4 Predictions for FAIR energies 



The CBM Collaboration at GSI is aiming at charmonium measurements at 
the future FAIR facility [73]. This opens up the possibility to explore the 
charmonium suppression mechanism at lower bombarding energies of about 
25 A-GeV in Au+Au collisions. First predictions for central reactions within 
the comover model have been reported in Ref. [22J. Here we extend the earlier 
studies to the full centrality dependence of the J/ty suppression also within 
the 'threshold scenario' and additionally provide predictions for the ^ to J/^> 
ratio. The corresponding HSD results are displayed in Fig. [10] for the survival 
probability Sj/q, (left plot) and ratio ^' to J/ty (right plot) as a function of the 
number of participants N part . The blue lines reflect the 'threshold scenario' for 
Ejjq, = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 , e*> = 6.55 GeV/fm 3 while the violet 
line stands for the 'threshold scenario' with a more realistic value of e%,i, i.e. 
Ejjq, = 16 GeV/fm 3 , e Xc = 2 GeV/fm 3 , e& = 2 GeV/fm 3 . The solid red lines 
denote the results for the comover absorption model with the standard matrix 
element squared |M | 2 = 0.18 fm 2 /GeV 2 . 

We note that in Au+Au reactions at 25 A-GeV the standard nuclear suppres- 
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Fig. 10. The survival probability Sj/^ (left plot) and ratio VP' to J/^ (right plot) as a 
function of the number of participants N par t i n Au+Au reactions at 25 A- GeV. The 
blue lines (with open dots) reflect the 'threshold scenario' for £j/$ = 16 GeV/fm 3 , 
e Xc = 2 GeV/fm 3 , e^/ = 6.55 GeV/fm 3 while the violet line (the lower line with 
open dots on the r.h.s.) stands for the 'threshold scenario' for Ej/q, = 16 GeV/fm 3 , 
e Xc = 2 GeV/fm 3 , e^i = 2 GeV/fm 3 . The solid red lines (full dots) denote the 
results for the comover absorption model with the standard matrix element squared 
|Mo| 2 = 0.18 fm 2 /GeV 2 . The dashed line (l.h.s.) represents the HSD calculations 
including only dissociation channels with nucleons. 

sion of J/ 1 ® (dashed line in the left plot) almost coincides with the 'threshold 
scenario' (solid line with open dots in the left plot) since only a very low am- 
mount of Xc and no J/^> are melted at the energy densities reached in these 
reactions. On the other hand the comover density decreases only moderately 
when stepping down in energy from 158 A- GeV to 25 A- GeV such that the 
J/ty survival probability in the comover absorption model (lower solid line in 
the left part) is substantially lower. This also holds for the to J/*f> ratio ver- 
sus centrality where even a lower threshold = 2 GeV/fm 3 leads to a ratio 
(middle line in the right part) that is clearly above the result achieved in the 
comover absorption model (lower line in the right part). Consequently the dif- 
ferent dissociation scenarios may well be distinguished in future charmonium 
measurements at FAIR. 



5 Summary 



In summarizing this work we have found that present data on charmonium 
suppression for Pb+Pb and In+In reactions at top SPS energies compare well 
with microscopic transport calculations in the comover model involving only 
a single parameter for the average matrix element squared | M 1 2 that fixes the 
strength of the charmonium cross sections with comovers. This holds for the 
J I ^ suppression as well as the ^' to J/ $ ratio versus collision centrality. The 



19 



bare 'QGP threshold scenario' gives satisfying results for the J/^> suppression 
for both systems at 158 A-GeV but fails in the VP' to J/^> ratio since too many 
ty' already melt away for a critical energy density of 2 GeV/fm 3 at 158 A-GeV. 
Only when assuming the to dissolve above ~ 6.5 GeV/fm 3 a reasonable 
description of all data is achieved in the 'QGP threshold scenario'; this thresh- 
old, however, is not in accordance with present lattice QCD calculations such 
that the 'threshold scenario' meets severe problems. This also holds for a com- 
bined model including 'threshold melting' and reduced 'comover absorption' 
as shown in Section 3. 

On the other hand the different scenarios can clearly be distinguished at FAIR 
energies (of about 25 A-GeV) where the centrality dependence of the J/\& sur- 
vival probability and the to J/ty ratio are significantly lower in the comover 
absorption model. This result comes about since the average comover density 
decreases only moderately with lower bombarding energy whereas the region 
in space-time with energy densities above critical values of e.g. 2 GeV/fm 3 
decreases rapidly and ceases to exist below about 20 A-GeV even in central 
collisions. This ideally might open up the possibility to measure excitation 
functions of the J/\& survival probability and the W to J/\& ratio in central 
Au+Au collisions, where clear steps would indicate the presence of 'melting 
thresholds' whereas a smooth excitation function would be in favor of the 
comover absorption approach. 
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